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Scanning tunneling microscopy and spectroscopy have been used to investigate the femtosecond 
dynamics of Dirac fermions in the topological insulator Bi 2 Se 3 ultrathin films. At two-dimensional 
limit, bulk electrons becomes quantized and the quantization can be controlled by film thickness at 
single quintuple layer level. By studying the spatial decay of standing waves (quasiparticle inter¬ 
ference patterns) off steps, we measure directly the energy and film thickness dependence of phase 
relaxation length 1^ and inelastic scattering lifetime r of topological surface-state electrons. We find 
that r exhibits a remarkable (E — Ef )~ 2 energy dependence and increases with film thickness. We 
show that the features revealed are typical for electron-electron scattering between surface and bulk 
states. 

PACS numbers: 68.37.Ef, 73.21.Fg, 73.50.Bk, 72.10.Fk 


Topological insulators (TIs) are generally character¬ 
ized by a bulk band gap and an odd number of topo¬ 
logically protected metallic surface states mu. Strong 
spin-orbit interaction brings about many exotic physi¬ 
cal properties in TIs and ensures the robustness of the 
surface metallic states against non-magnetic impurities. 
With the nontrivial spin textures of surface states, three- 
dimensional (3D) TIs in reduced dimensions have gen¬ 
erated tremendous interests since their discovery, due 
to the enhanced surface-to-volume ratio and potential 
for developing future dissipationless spintronics USO]. 
Angle-resolved photoemission spectroscopy (ARPES) on 
Bi 2 Se 3 films has revealed an energy gap opening associ¬ 
ated with the hybridization of top and bottom surface 
states in the two dimensional (2D) limit [5]. Extensive 
transport experiments have demonstrated thickness de¬ 
pendent electrical resistivity and magnetoresistance in 
ultrathin films of TIs , which has been suggested as the 
combined effects of strong electron-electron (e-e) interac¬ 
tions and topological delocalization (SHU- On the other 
hand, ARPES can act as a powerful tool to probe directly 
the dynamics of Dirac fermions m • So far as we know, 
only a few such experiments were reported, but focused 
primarily on electron-phonon (e-p) and electron-impurity 
scattering and led to inconsistent results EHH]. Thus 
far, the e-e interactions of topological surface states in 
TIs have scarcely been addressed separately and estab¬ 
lished conclusively. 

The phase relaxation length Z^, defined as the average 
distance a quasiparticle can propagate without losing its 
phase memory, is usually thought of as a key quantity to 
describe the femtosecond dynamics of electrons in solid 
state physics including TIs. Therefore an accurate esti¬ 
mate of 1(f) in TIs is highly desirable, but challenging in 
experiment. In B^Ses films l$ at the Fermi level (Ep) 
was recently studied by magnetoresistance measurements 


IMS], but the obtained results of the intrinsic l^Ep) as 
well as its thickness dependence are not consistent. This 
may stem from varying surface defect scattering (i.e. step 
edges and residual impurities) in the samples from differ¬ 
ent groups. An alternative access to inelastic lifetime r, 
and hence to /</> = vpr (vp is the group velocity of elec¬ 
trons) , has recently become available by analyzing quan¬ 
titatively the energy widths of the quantized topologi¬ 
cal surface states with a scanning tunneling microscope 
(stm) ng. Such local technique eliminates the residual 
scattering and allows the estimate of 1$ at various ener¬ 
gies. However, the method demands elaborate modeling 
to remove the instrumental effects, which also severely 
limits an accurate determination of 1$. 

In this Letter, we study the decay of standing waves 
(SWs) off descending straight step edges in B^Ses films 
[Fig. 1(a)] with scanning tunneling microscopy and spec¬ 
troscopy (STS). We measure directly l$ (or equivalently 
r) as a function of energy and film thickness. All ex¬ 
periments are carried out on a Unisoku ultra-high vac¬ 
uum STM system equipped with molecular beam epitaxy 
(MBE) for film growth. High-quality B^Ses films with 
controlled thickness are grown by co-evaporating high- 
purity Bi (99.999%) and Se (99.999%) onto graphitized 
677-SiC(OOOl) substrate kept at 220°C [5j H6j . except 
for 19 quintuple layer (QL) film at 260°C. Prior to the 
STM/STS measurements at 4.8 K, a polycrystalline W 
tip is cleaned by electron-beam heating in MBE cham¬ 
ber, and then calibrated with STS of Pb films grown on 
a Si(lll) substrate in an energy range of -1.0 eV ^ 1.5 
eV [T7] . All the differential conductance dl/dV spectra 
and maps are acquired using a standard lock-in technique 
with a small bias modulation of 10 mV at 987.5 Hz. 

Figure 1(b) depicts a series of differential conductance 
dl/dV spectra in the as-grown B^Ses films with various 
thicknesses ranging from 1 QL up to 50 QL. Above ~10 
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FIG. 1. (a) STM topographic image of Bi 2 Se 3 film with a 
nominal thickness of 15 QL ( V s — 4.0 V, I = 0.1 nA, 200 
nm x 200 nm). The step edges are invariably oriented along 
the closed-packed directions of Bi 2 Se 3 , indicated by the three 
red dashes, (b) A series of dl/dV curves showing the quan¬ 
tum quantization of electrons (triangles) in the 2D limit. The 
spectra have been vertically offset for clarity. Set point: V s 
= 0.2 V, I = 0.1 nA. (c) Inverse of the energy separation be¬ 
tween the quantized states marked by red and blue triangles 
in (b) as a function of film thickness t. 

QL, the spectra display a minimum at ^0.2 eV, which 
corresponds to the Dirac point of topological surface 
states and agrees with those reported previously uses. 
As the films reach the 2D limit, however, additional dis¬ 
crete states (marked by colored triangles) are discernible 
in the valence band and alter in energy position with film 
thickness t. These are reminiscent of the M-shaped 2D 
electron gas observed by ARPES Gang, which has been 
contentiously interpreted as quantum well states (QWSs) 
0 HU HQ or the expansion of van der Waals’ spacing 
caused by impurity intercalation [22] . Our experiments 
are conducted on MBE-grown clean films and little suf¬ 
fers from impurity intercalation, and the results supports 
strongly the QWSs scenario. This is further evidenced by 
revealing the film thickness t dependent energy separa¬ 
tion (A) between adjacent QWSs, plotted in Fig. 1(c): 
the linear relationship between 1/A and t demonstrates 
indubitably that all the discrete states observed here and 
also by previous ARPES experiments originate merely 
from the quantization of bulk electron states [23 . Con¬ 
sidering the vanishing QWSs in thick B^Ses films [Fig. 
1(b)], we note that the quantization of electron in B^Ses 
films and bulk single crystals are driven distinctively and 
primarily by the dimensionality and band-bending effects 
HU HU El], respectively. 

In order to measure the phase coherence length we 
select straight step edges in an almost defect-free region 
(^100 nm x 100 nm) to minimize unwanted electron 
scattering. On the upper terrace, we acquire simulta¬ 
neously the topographic image (shown in Fig. SI) [24] 
and differential conductance dl/dV map in an area of 
30 nm x 30 nm, similar to the previous reports [25] [26] . 
Here care should be taken to determine the sample lo¬ 
cal density of states (LDOS) from which is extracted. 
In the vicinity of step edges, strong spatial oscillation of 
SWs alters the tip-sample distance, and consequently the 
measured dl/dV will not exactly represent the LDOS. 


To estimate more accurately, we have recovered the 
LDOS from STM topography and dl/dV map using a 
well-tested algorithm m, as illustrated in Figs. 2(a-f) 
and 3(a-e). Here the work functions of W tip <p t = 4.55 
eV and B^Ses <j) s = 5.5 eV [28] are adopted. Further 
variation of (f) t + (j) s up to 0.5 eV is found to little af¬ 
fect the reconstructed LDOS. Moreover, a significantly 
small lock-in bias modulation of 10 meV was adopted to 
reduce its additional decay in dl/dV map, which may 
complicate the estimate of [29]. By slightly rotating 
the LDOS maps to align the step edge vertically, we av¬ 
erage the LDOS as a function of distance x from the step 
edge, as illustrated in Figs. 2(g-l) and 3(f-j). 

As expected, the SWs are clearly visible on the terrace 
and decay as a function of x. In ordinary 2D electron gas 
[29l [3T] , it has already been demonstrated that the decay 
is primarily governed by the combined effects of geomet¬ 
ric factor (~l/x 0 - 5 ) and phase decoherence 
whereas the latter has completely been neglected in TIs 
by assuming 1$ = oc mm- Our attempts to fit the 
LDOS by the predicted cos(qx + (j))/x a fail if the loss of 
phase coherence is not considered (red curves in Figs. 2 
and 3), in particular for those at energies far away from 
Ep. The observed LDOS invariably exhibit a faster de¬ 
cay than l/x a , where a = 0.5 and 1.5 are found to best 
represent the LDOS for t = 1 and for t > 5, respec¬ 
tively (here q and (j) denote the scattering wave vector 
and phase shift, respectively). Our careful examination 
shows better fits of the LDOS if the phase decoherence 
effect is included (blue curves in Figs. 2 and 3), namely 

LDOS(q, x) oc |r| cos (^ + $ e ~*/4, (1) 

where r indicates the electron reflection coefficient by 
the step-edge barrier m- For a finite barrier (e.g. step 
edge), \r\ is generally smaller than 1, but it does not affect 
the estimate of q and 1$. Our observation thus presents 
a compelling evidence that the phase coherence has been 
appreciably lost in the electron scattering process. 

Now we comment on the different geometric param¬ 
eter a between 1 QL and thicker (t ^ 5) Bi 2 Se 3 films. 
In TIs with weak warping effect, such as B^Ses investi¬ 
gated here, the unusual topological spin texture prohibits 
backscattering, leading to a faster decay of SWs off steps 
as 1/x 1-5 [26] 03], in stark contrast to the 1/x 0,5 decay 
in an ordinary 2D electron gas El HD]- This matches 
excellently with our observations for thicker B^Ses films 
in Figs. 3 (a-j). For single QL B^Ses film, however, the 
top and bottom surface states can talk with each other 
and break the momentum-spin helical locking symme¬ 
try of topological surface states, confirmed recently by 
spin-polarized ARPES experiments [341 35]. As a con¬ 
sequence, the backscattering from the top surface state 
is allowed, leading to the 1/x 0,5 decay in single QL 
Bi 2 Se 3 film [Fig. 2], similar to the ordinary 2D electron 

















3 




r 

m 


(f) 


x (nm) 


x (nm) 


x (nm) 


x (nm) 


x (nm) 


x (nm) 


FIG. 2. (a-f) LDOS maps of 1 QL Bi 2 Se 3 film taken on the upper terrace near a step edge (I = 0.1 nA, 30 nm x 30 nm), and 
(g-1) the averaged LDOS (black dots) as a function of distance x from the step edge for various energies. Note that the mean 
value of each LDOS map has been subtracted. The SWs patterns with different decay rate are clearly evident. Blue and red 
lines show the best fits of the decaying LDOS with and without a loss of phase coherence considered, respectively. 
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FIG. 3. (a-e) LDOS maps (/ = 0.1 nA, 30 nm x 30 nm), and (f-j) the averaged LDOS (black dots) as a function of distance x 
from the step edge for various film thicknesses and energies. Blue and red lines show the best fits of LDOS with and without 
a loss of phase coherence considered, respectively, (k) Energy E- q dispersion relation deduced by fitting the averaged LDOS 
around the step edges, with q indicating the scattering wave vector. The Dirac points for Bi 2 Se 3 films from 1 QL to 200 QL 
are estimated to be -0.51 eV, -0.31 eV, -0.30 eV, -0.16 eV and -0.08 eV by intercept of the E- q dispersions with the energy axis 
at q = 0. 


gas. Therefore, our simple scattering experiments off 
steps have unequivocally demonstrated not only that the 
backscattering of topological surface states by nonmag¬ 
netic impurities (e.g. step edges) is indeed suppressed 
[25} 36, but also that the top and bottom surface states 
hybridize strongly and break the chiral spin texture in 




FIG. 4. (a) Inelastic lifetime r of topological Dirac fermions 
in Bi 2 Se 3 as a function of energy and film thickness. Colored 
lines depict the best fits of lifetime to r = A (E — Ef )~ 2 for 
various film thicknesses, (b) Double-logarithmic plot of the 
extracted A from (a) versus film thickness t. A increases as 
£0.4i±0.08 thickness t. 


the 2D limit [5j [34} [35] . More significantly, the phase de¬ 
coherence effect on the SWs decay is well separated from 
the geometric factor in Eq. (1), allowing a more accurate 
absolute estimate of 1$. Figure 3(k) summarizes the E- 
q dispersion relation of B^Ses surface states, where the 
wave vector q are extracted from the fits of LDOS to Eq. 
(1). Colored lines reveal the linear fits to the data for 
various film thicknesses. The approximately equal slope 
hints at the same dominant scattering channel for vari¬ 
ous film thicknesses. Using q = 2k [26] . we deduce the 
Fermi velocity to be vp = 4.0 =b 0.2 eV-A along the Y-M 
direction. 

Figure 4(a) plots the inelastic lifetime r, calculated 
via r = 1$/vp, as a function of energy and film thickness. 
The error bars primarily indicate the standard derivation 
of r values obtained from the LDOS fits over various x 
ranges. Note that in our measurements we have used a 
typically small tunneling current of I = 0.1 nA, at which 
one excited electron about every 1.6 ns (^> r) is injected 
into (V s > 0) or knocked out of (V s < 0) the samples. 
This avoids nonequilibrium probe and electron cascade 
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effect. It is clearly evident from Fig. 4(a) that the inelas¬ 
tic lifetimes r of topological Dirac fermions are strongly 
reduced when the energy E departs away from Ep- This 
is typical for e-e scattering, since the e-p lifetimes are 
generally little dependent of quasiparticle energy for the 
small energy range studied here. Moreover, we find that 
r changes with energy in a power-law manner. In Fig. 
4(a) the five colored lines shows the best fits of r for 
various B^Ses films by 


Such power law behavior agrees well with 3D interacting 
Fermi liquid theory for electron-hole pair creation [37] , 
similar to previous observation for conventional surface 
states on Ag(lll) and Cu(lll) [29] . We thus believe 
that the phase relaxation lengths l$ or lifetimes r mea¬ 
sured here are predominantly governed by inelastic e-e 
scattering in B^Ses films. This is understandable since 
e-p scattering is expected to be insignificant at the low 
temperature of 4.8 K. 

Finally we shed more insight into the observed e-e scat¬ 
tering by studying the film thickness t dependence of A, 
plotted in Fig. 4(b). Here a larger A means a longer phase 
coherence length or lifetime r. It is found that A in¬ 
creases with film thickness as A oc £°- 41 ± 0 - 08 . if the inelas¬ 
tic scattering proceeds mostly via the 2D channel, i.e. the 
intraband scattering within the surface states, A would 
depend little on t since the surface should not alter with 
film thickness. Therefore, the film thickness-dependent 
inelastic lifetime r (or A) reveals directly a predominant 
3D scattering channel in B^Ses TIs, namely e-e scatter¬ 
ing between surface and bulk states. This is consistent 
with the power law dependence of r on energy observed 
above, predicted in 3D interacting Fermi liquid theory. 
The enhanced e-e scattering in thinner film is primarily 
related with reduced dimensionality, which results into 
weaker electron screening and hence stronger Coulomb 
interaction between electrons. As compared to the e- 
e scattering in Ag and Cu m , A appears substantially 
smaller in B^Ses, indicative of weaker electron screening 
in TIs than that in ordinary metals. 

In summary, our detailed STM and STS study of the 
SWs decay in MBE-grown B^Ses films has revealed the 
significant loss of phase coherence. By fitting the LDOS, 
we have estimated the energy and film thickness depen¬ 
dence of the phase coherence length or inelastic lifetime 
r. The inelastic lifetime r depends on both energy and 
film thickness in a power law way. This has convincingly 
demonstrated the significance of e-e scattering between 
surface and bulk states in TIs. The present STM/STS 
approach avoids macroscopic integral measurements, al¬ 
lowing for a better estimate of and r. It should be 
highly desirable to generalize this technique to other TIs, 
and also to probe the temperature dependence of SWs 


from which the contributions of e-p scattering on the 
phase decoherence can be quantified. 
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